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Abstract 

The one-dimensional totally asymmetric simple exclusion process (TASEP) is 
considered. We study the time evolution property of a tagged particle in the TASEP 
with the step initial condition. Calculated is the multi-time joint distribution function 
of its position. Using the relation of the dynamics of the TASEP to the Schur process, 
we show that the function is represented as the Fredholm determinant. We also study 
the scaling limit. The universality of the largest eigenvalue in the random matrix 
theory is realized in the limit. When the hopping rates of all particles are the same, 
it is found that the joint distribution function converges to that of the Airy process 
after the time at which the particle begins to move. On the other hand, when there 
are several particles with small hopping rate in front of a tagged particle, the limiting 
process changes at a certain time from the Airy process to the process of the largest 
eigenvalue in the Hermitian multi-matrix model with external sources. 
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1 Introduction 



Dynamics of a nonequilibrium system is one of the most stimulating topics in statistical 
mechanics. The dynamical property is determined by the interplay among various elements 
such as interaction, initial and boundary conditions, and so on. But little is known about 
how it depends on these elements and realizes universality in a situation where the elements 
are intertwined. The aim of this article is to approach such questions by analyzing the 
dynamics of the one-dimensional asymmetric simple exclusion process (ASEP). 

The ASEP is one of the most typical models of interacting particle processes where 
particles diffuse to a preferred direction with hard core repulsive interaction [U |2[ |3l IH |5] . 
In spite of the simplicity of the model, it has been known that the ASEP shows various 
interesting phenomena caused by the collaboration between the diffusion and exclusion 
effect. Furthermore, in the one-dimensional case, it has an integrable mathematical struc- 
ture which allows us to analyze some physical quantities exactly. For example, in the 
open boundary condition a steady state can be constructed using the matrix product 
method and g-orthogonal polynomials, and the boundary-induced phase transition can be 
discussed [6l[71[8]. 

On the other hand, the dynamical (non-stationary) properties of the one-dimensional 
ASEP have been also analyzed recently. Among various studies on the topic, we focus 
on the problem of the diffusion of a particular particle in the one-dimensional ASEP. (We 
call the particle a tagged particle.) This is a fundamental and elementary problem for the 
understanding of the dynamics of the ASEP. We are interested in how initial conditions 
and the exclusion effect through hard-core repulsive interaction affect the diffusion property 
and how it is different from the (normal) diffusion of the Brownian particle. 

The scaling exponent of the diffusion provides insights into these problems. It was 
found that it depends on the initial conditions. In a steady state with a given density, 
the position fluctuation of a tagged particle grows as t^/^ [9]. The exponent 1/2 is the 
same as that of the Brownian particle. On the other hand, for the fixed initial conditions 
where initially we create particles with a given density and fix the initial configuration, 
the exponent changes from 1/2 to 1/3 pUl [TT] . Such a diffusion is called an anomalous 
diffusion. This exponent 1/3 can be understood from the fact that the time evolution of 
the one-dimensional ASEP can be described by the one-dimensional Kardar-Parisi- Zhang 
universality class p!2] . In the study of the current in the one-dimensional ASEP on a ring, 
the exponent 1/3 was calculated by the Bethe ansatz technique [131 HI]. 

Note that the two cases above have a common property that the density is invariant 
under both time and space translation and a tagged particle moves with a fixed velocity on 
average. On the other hand, let us now consider another initial condition where all sites on 
the left of a certain bond are occupied while all right sites are empty. We call this initial 
condition a step initial condition. This is the typical initial condition where the density 
and the average velocity of a tagged particle change with time. Under the initial condition, 
we focus on the Mth particle from the right as a tagged particle and consider the time 
evolution of the particle when M is large (but finite). In order to study this problem, 
we investigate the limiting behavior of the diffusion of the tagged particle as both M and 
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time go to infinity. In the totally asymmetric simple exclusion process (TASEP), where 
particles move only to the right, the average position of a tagged particle has been given 
in the scaling limit [15]. On the other hand, we are interested in the limiting process itself 
which the tagged particle obeys, as well as the average position and the scaling exponent. 

Such an attempt as getting more sophisticated information than the scaling exponent 
has been recently discussed in the study on the current fluctuations in the TASEP. Not 
only the exponent but also the scaling function were obtained and it was revealed that the 
function is equivalent to the largest eigenvalue distribution in the random matrix theory 
with various universality classes. This development is based on the results on the longest 
increasing subsequence in a random permutation [Ml |l7l [181 UHl EO] . 

Remarkable is that the scaling function detects a difference in initial conditions even 
if the scaling exponent does not. For instance, in the case of the step initial condition, 
the scaling function is equivalent to the GUE Tracy- Widom distribution [21], which is the 
largest eigenvalue distribution of the Gaussian unitary ensemble (GUE) [TTl [221 ES]- In 
other initial conditions, on the other hand, these are described as the largest eigenvalue 
distributions with other universality classes in RMT such as Gaussian orthogonal ensemble 
(GOE) [21], two independent GOE's which is denoted by GOE^ [191 [20] , and so on [221 [25] . 
In addition, the equal-time multipoint distribution of the current fluctuations is obtained in 
the step and periodic initial conditions ^26j by use of the Green function of the TASEP [27] . 

In this article, we calculate the multi-time distribution function of position fluctuations 
of a tagged particle in the TASEP with the step initial condition based on the techniques 
of the directed polymer problem in a 01 matrix [281 ISH [301 EI] ^-nd stochastic growth of 
Young diagram characterized by the Schur process [321 [331 EH [351 ESI [37]. We express 
the function as the Fredholm determinant and discuss the scaling limit by the saddle point 
analysis. We remark that what will be considered in this paper is the pointwise convergence 
of the kernel. The convergence of the Fredholm determinant itself is certainly expected 
to hold from previous works on related problems [381 139] , but more elaborate asymptotic 
analysis is necessary for its proof. In the study on the ASEP, such correlations between 
different times have not been studied yet although the Green function [231 l2ll SO] and the 
equal-time joint distributions of particle positions [SHI E] have been studied recently. In 
the one-time case, on the other hand, the position fluctuation is essentially the same as the 
current fluctuation which has already been studied in [TTl [25] . But it has not been stated 
as the position fluctuation. 

Furthermore, we also consider the dependence of the distribution function on defect 
particles located in front of the tagged particle. Here defect particle means a particle with 
different hopping rate from other normal particles and we assume the number of them is 
finite. Such a situation as the TASEP with particle-dependent hopping rates has been 
recently discussed in [^ |33l Hlj . We clarify that when their hopping rate is smaller than 
that of normal particles, the limiting behavior of the multi-time distribution changes due 
to their presence although their number is finite. 

The paper is organized as follows. In the next section, we describe the definition of the 
model. The main result of this article is given in Section [31 In Section [H we express the 
multi-time distribution of position fluctuations in the TASEP as Fredholm determinant. 
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For this purpose, we first show that the directed polymer problem on a 01 matrix is related 
to the time evolution of a tagged particle in the TASEP with the step initial condition. 
Next, by applying the dual Robinson-Schensted-Knuth algorithm to the 01 matrix, we 
map the tagged particle problem in the TASEP to the stochastic growth of Young diagram 
described by the Schur process. In Section [5], we discuss the two types of scaling limit 
for the multi-time distribution function, the case where both time t and the label of the 
tagged particle M tend to infinity, and the one where t goes to infinity with M fixed. Some 
discussions and concluding remark are given in Section [6] and Section [7] respectively. 

2 Model 

In this section, we define the model and the quantity which we study in this article. Let us 
consider the one-dimensional infinite lattice and particles as illustrated in Fig. 1. Each site 
can be occupied by at most one particle. Suppose all particles are labeled an integer i from 
the right. A particle moves stochastically obeying the following rules. During each time 
step between t G {0, 1, ■ ■ ■ } and t + 1, the particle labeled i hops to the right neighboring 
site with probability 1 — gj (0 < < 1) and stays at the same site with probability if the 
right neighboring site is empty. On the other hand, if the site is occupied by the particle 
labeled i — the ith particle stays at the same site with probability 1. This incorporates the 
exclusion effect which describes the hard core repulsive interaction. When hopping rates 
of all particles are the same, i.e., g^ = g(Vi), the system is the usual totally asymmetric 
simple exclusion process (TASEP) with parallel update. We allow the particle-dependent 
hopping rates and study the effects of defect particles. 

In this article, we consider the step initial condition in which all sites on the left of a 
certain bond between two sites are occupied by particles and all sites on the right of the 
bond are empty. Fig. 1(b) illustrates the initial condition. 

Under these settings, we label the particles at 1,2, ■■ ■ from the right (see Fig. 1(b)) 
and consider the dynamics of the particle labeled M. (This particle which we focus on is 
called the tagged particle.) For this purpose, we set the position coordinate such that the 
tagged particle is at the origin at time as depicted in Fig. 1(b). Let us define L{t, M) as 
the position of the Mth particle at time t. In other words, L(t, M) represents the distance 
travelled by the Mth particle from time through t. 

The quantity which we will investigate is the multi-time joint distribution function for 
L{t, M) defined as 

Prob (L(ti, M) > £i, L(t2, M) > £2, ■ ■ ■ , ^(^m, M) > £„) . (2.1) 

In particular, we are interested in the asymptotic behavior of (12. ip as t goes to infinity. In 
this article, two types of scaling limit will be discussed. The first one is the limit where 
both time t and the label of the tagged particle M go to infinity with the ratio t/M fixed. 
In the second one, on the other hand, we take the t —>■ 00 limit with M fixed. 

In order to see how the whole particles from the first to Mth travel, we performed 
the numerical simulation of which the result is shown in Fig. 2. We considered the time 
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evolution of the particles from the first to 100th and assumed the hopping rate of all 
particles is 0.9 except the first, 25th, 50th and 75th particles whose hopping rate is 0.8. 
Thus these four particles represent the defect particles. The typical example of the time 
evolution of the 100 particles from time through 3000 is illustrated in Fig. 2(a). In this 
figure, the vertical axis represents the position coordinate introduced in Fig. 1(b), whereas 
the horizontal axis is time axis, and each particle is expressed as +. From the macroscopic 
point of view such as Fig. 2(a), we only find that they travel in a group. However, if we 
zoom the configuration of the particles in a certain time zone, we can see the microscopic 
pattern of the configuration. Here we show two characteristic configurations which are 
depicted in Fig. 2(b) and 2(c). Fig. 2(b) (resp.(c)) represents the configuration of the 
particles when time is around 200 (resp. 3000). When the time zone which we focus on is 
early enough as in Fig. 2(b), the particles move forming only one group and the four defect 
particles are included in the group. On the other hand, when time has passed sufficiently 
as in Fig. 2(c), they form the four groups automatically and the top of each group is the 
defect particle. The similar platoon structure was also discussed in the study of the TASEP 
with disordered hopping rates [l5l HHl H?] . In fact, as will be explained in the next section, 
the position fluctuation of a tagged particle also changes with time corresponding to the 
change of the configuration. 



3 Main results 

3.1 Multi-time distributions 

Our whole discussions in this article are based on a fact that the tagged particle problem 
for the step initial condition has a nice combinatorial structure which admits us to obtain 
a closed expression for the multi-time joint distribution in the form of the Fredholm deter- 
minant. The results are summarized in the following theorem. The proof will be presented 
in Section m 

Note that due to the rule of the TASEP, the particle labeled M cannot move when 
t < M. Thus, throughout the article, we consider the multi-time function (12.11) under the 
condition 

ti > M. (3.1) 



Theorem 1 

When dSH) holds, 

Prob {L{t^, M) > £i, L(t2, M) > 4, ■ ■ ■ , L{t^. M) > C) = det (1 + Kg) . (3.2) 
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Here det(l + Kg) is the Fredholm determinant defined as 
det (1 + Kg) 

oo ^ m oo m oo 

= X1^5Z Yl '"Yl 9{W,xi)---g{tn,;Xk)det{K{tni,xi;tn,nXi'))li,^^, 

k=0 711=1 xi=—oo nf^=lxk=—oo 

(3.3) 

where 

g{tn; x) = -X{t,,-M+i-e,,,oo)ix), (ra = 1, 2, . . . , m), 

«"'"W = (n '{r"'^'- (3.4) 
I 0, otherwise, 

and the kernel K{ti,xi;t2,X2) is given by 

K{ti,xi;t2,X2) = K{ti,xi;t2,X2) - (f)t,,t2ixi,X2), (3.5) 

k(ti,Xi;t2,X2) 

1 f dzi [ dz2 z, (1 + 1/^2)^^-^+^ TT ^ ~ ~ ^^-"^ -"2^ ^6^ 
(27rz)2 y^^^ ^1 Jc^^ Z2 z, -^2(1 + l/^i)*-^^+i Al 1 _ _ g^;,^ ' ^ • ^ 

4>t^,t2{xi,X2) = \^^''^^ ^ ^ ^' (3.7) 

I U, t\ >_ 12- 



Here Cr in f l3.6p and (13.71) denotes a contour enclosing the origin anticlockwise with radius 
R and Ri {i = 1, 2) in (13. 6p satisfy the conditions R2 < Ri and 1 < fZi < (1 — gj)/gj. 



3.2 Scaling limit 1 (M ^ 00) 

Using the results in Theorem 1, one can study the asymptotics of the joint distribution. 
In this subsection, we consider the scahng hmit such that both t and M go to infinity with 
their ratio 

u = t/M (3.8) 

fixed, u represents the scaled time. Notice that we also take the limit of the label of the 
tagged particle M. This may sound strange if the particle we are focusing on is varied as M 
goes to infinity. Rather we think that the results below give the asymptotic behaviors of the 
path statistics of the tagged particle for large but finite M which is fixed. Alternatively one 
might be able to interpret them as giving the asymptotic behaviors of the path statistics 
of the tagged particle near a large but finite t. 

We also discuss the effect of defect particles. In particular, we consider the situation 
where there are finite n defect particles with stay rates qi{i = 1, ■ ■ ■ ,n) in front of the 
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tagged particle. For this purpose, we set the stay rate qi of ith. particle as follows. For 
some set {aj}i=i,... ,n C {1,2, ■ ■ -M} with n fixed, we assume 

if i = a,-, , ^ 

^' (3.9) 
otherwise. 

Note that the particle with the rate qj represents the defect particle whose label is aj 
whereas the remaining particles with the rate q are the normal particles. 

3.2.1 Average position 

The average position of the tagged particle divided by M has a deterministic limit as 
M ^ oo. Let us call it A{u), 

A{u) = lim ^ TT^- (3-10) 

Al ^oo IVl 

We also define q = maxjgj}. (Note that q is the stay rate of the slowest defect particle.) 
If g > g, which is the case where the slowest defect particle is slower than the normal 
particles, we obtain 

A{u)= <A2{u), if Y^<M<M„ (3.11) 
[Ag{u), if Uc < u, 

where Uc = {(f — 2qq + q)/{q — qY and 

A2{u) = (1 - q)u - (1 - 2q) - 2^q{l - q){u - 1), (3.12) 

AG{u) = {l~q)u-^-^^^. (3.13) 

q-q 

A2{u) can be understood from Theorem 1.1. in [17J. If g < g, the average position is 
represented as only the first and second cases in (13. lip , {uc and Ag{u) do not appear in 
the case.) Note that A2{u) does not depend on the stay rate g of the slowest defect particle 
whereas Ag{u) does. A{u) is illustrated in Fig. 3. 

Next we consider the position fluctuations of the tagged particle around the average 
position. There are four typical regions shown in Fig. 3. In each region, we can obtain 
the result on the position fluctuations by focusing on a point within the region and taking 
a proper scaling around the point. Note that in Fig. 3, the regions 2 and 4 spread in 
the form of a line whereas the regions 1 and 3 are point-like. Thus, in the region l(resp. 
3), we always consider the fluctuation property around the point 1/(1 — g) (resp. Uc = 
if -2qq + q)/{q- qf). 
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3.2.2 Region 1 (m = 1/(1 - g)) 

This is the region where the tagged particle just starts to move. This corresponds to 
u = 1/(1 - q). We set 

M 1 

ti = - + DiM2T,, (3.14) 

1-q 

where Di = y/q/{l — q). Note that the time ti (13.141) is approximately equal to the time 
on which some leftmost holes arrive at the site occupied by the tagged particle. One easily 
finds that the arrival time can be scaled as in (13.141) which is the same scaling as the central 
limit theorem. 

Applying the scaling to Theorem 1, we have the theorem as follows. The proof will be 
given in Section [5?T1 

Theorem 2-1 

lim Prob (L(ti, M) > 4, L{t2, M) > £2, ■ " " , ^(^m, M) > £„) = det (1 + /Ci^i) . (3.15) 

M— >oo 

Here the Fredholm determinant in the right hand side is defined in (13.31) where gi{Ti,x) = 
— X(_oo/i)(a;). X(a,b){x) is defined in (13.41) . The kernel is given by 

^ . X _ /Em=oV^i(a^i- ^1)^2(^2 -m,r2), n>T2, 

[-2^m=0'^n^^ +fn+l,Ti)lp2{X2+fn + l,T2), Ti < T2. 

Here iIji{xi,Ti) can be represented by use of the parabolic cylinder function Dn{x) in fj^ 
as 

^i(xi,ri) = — / dze--^^'z^^'' = -j=D^,.,{t), (3.17) 

^TTi y_joo+e V27r 

where e is taken to be positive, and 

U^.r2) ^ ^ I 4^e-^-- ^ I "^^^ (^) ' ^ °' (3.18) 



2tiiJc,w-^+^ [o, X2<0 

where Ci in (I3.18P represents the contour enclosing the origin anticlockwise with radius 1, 
and Hn{x) is the Hermite polynomial with degree n fj^ . 

The special case where ti = ■ ■ ■ = = has appeared in the distribution of the height 
fluctuation in the "critical regime" of the oriented digital boiling model [30]. Recently the 
distribution in one time case (m = 1) has been also given in [49j. Our formula above gives 
a generalization to the multi-time version. 

In particular, the limiting distribution in the case m = 1 can be described as 

hm Prob {L{t, M)>i) = Y, E " " " E ^' ^''))n'=i , (3-19) 

k={) xi=0 a;fe=0 



8 



where 

K:i{t,Xi;t,X2) = - m,T)'ilj2ix2 -m,T). (3.20) 

m=0 

Note that in fl3.19l) the summation on k where k > i + 1 can be omitted. Thus in this 
case, we can obtain exactly the hmiting probabihty by calculating the determinants with 
finite rank. Another expression of the kernel (13.201) has recently appeared as the "discrete 
Hermite kernel" in [l9]. In [30j, the specific values of the probability are given in the case 
r = 0. 

3.2.3 Region 2 (1/(1 - q) < u < u^) 

This is the region where the effect of the defect particles does not affect the time evolution 
of the tagged particle and the dynamics of the ordinary TASEP is dominant. The typical 
picture of the time evolution of the whole particles from 1st to Mth in this region is 
illustrated in Fig. 2(b). In this figure, we are interested in the position fiuctuations of the 
bottommost particle. 
Let us scale as 

tj = uM + C{u)Mhj, (3.21) 
ij = A2{uj)M - D{u)Mhi, (3.22) 

where 1/(1 — g) < u < Uc and Uj = tj/M. A2{u) is defined in (13.121) and 

1 

^(") - 2(« - 1)« (l + ' (V^^ - ^ ' . (3.23) 

2 

D(u) = (. - 1)J,4(1 - ,)i H- ^/III^ ' - . (3.24) 

Note that the scaling exponent 1/3 in fl3.22p also appears in an anomalous diffusion of 
the tagged particle [TOldT]. The scaling exponents 2/3 and 1/3 in (13.211) and (I3.22p are 
characteristic of the one-dimensional KPZ universality class. They are expected to be 
universal for all models belonging to the KPZ universality class. 

Under the scaling defined above, we get the following theorem for the scaling function. 
This function is also universal but may depend on various elements such as initial conditions 
and the effect of the defect particles. For example, as we will show in Theorem 2-3, the 
scaling function in the region 3 is different from that in region 2 although we take the 
same scaling (13.211) and (13.221) in both regions. Thus, the function provides more detailed 
viewpoint in the KPZ universality than the scaling exponents. The proof of this theorem 
will be given in Section 15. 2[ 

Theorem 2-2 

lim Prob (L(ti, M) > ^i, L(t2, M) > 4, ■ ■ ■ , ^(^m, M) > £„) = det (1 + /Ca^) . (3.25) 
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Here the right hand side is the Fredholm determinant defined as 
det (1 + /Ca^) 

CO ^ m /•cxD ^ POO 

fc=0 ni=W-°° nfe=l-^-°° 

(3.26) 

where Qijj,^) (j = 1, ■ ■ ■ ,m) are defined in terms of X(a,b){x) fl3.4p . as 

6^(^j,0 = -X(s„oo)(0 (j = l,---,H, (3.27) 
anc? the kernel }C2{ti,C^i]T2, (,2) is given by 

r(rf-^f\ ,/„°°<iAe-*'^'""'AiK, + A)Ai(& + A), n > t^. 



The kernel /C2 is called the extended Airy kernel pO\ [51] . The process characterized 
by the Fredholm determinant with this kernel is called the Airy process |3H1 [52]. This 
process appears as the limiting process of the largest eigenvalue in Dyson's Brownian 
motion model [53] of the unitary class. The model is described as N x N Hermitian 
matrix where each independent element of H obeys the Ornstein-Uhlenbeck process. The 
transition probability density PiHf, Hj; r) from matrix Hi to Hj during r is given by 

P {Hf, H,; t) = Z.exp ^f^~!,7^^^' j , (3.29) 

where Zr is the normalization constant. If we choose the initial matrix Ho to be GUE 
random matrix, the joint density function of the probability that matrix Hj is at time tj 
is represented as the Hermitian multi-matrix model, 

P {Hi,ti] ■ ■ ■ ; Hm, tm) 



where Z is the normalization constant. Let Z*^*-* be the largest eigenvalue of Hi and we 
consider the quantity, 

Prob (/(^) < ai, • ■ ■ , < am) ■ (3.31) 

When we set 

U = ai = V2N + — (3.32) 

the limiting distribution (13.311) is described as 

lim Prob (/(I) < ai, ■ ■ ■ , l^""^ < am) = det (1 + /Ca^) . (3.33) 
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In the case m = 1, the distribution is the GUE Tracy- Widom distribution [2T] . 

Hence our Theorem 2-2 says that in the appropriate scahng hmit the dynamics of the 
tagged particle in this region is equivalent to the dynamics of the largest eigenvalue of 
Dyson's Brownian motion of the unitary type. In the context of the one dimensional KPZ 
universality class, the Airy process has already appeared in the study of the polynuclear 
growth (PNG) model. In [381 [52], the process has first appeared as the equal-time multi- 
point height fluctuation. Recently in [31], it has been shown that the process also describes 
the correlations among any "space-like" points in the space-time plane. Our result is closely 
related to the latter situation. 

3.2.4 Region 3 {u = Uc) 

This region is the border between the region 2 where the normal particles in the TASEP are 
dominant, and the region 4 where only the finite number of defect particles are dominant. 

First we consider the case where the stay rates qj {j = 1,2, ■ ■ ■ , n) of the defect particles 
are distinct. Taking the same scaling as ( I3.2ip - (13.24p with u = Uc, we obtain the following 
theorem. 

Theorem 2-3 



lim Prob M) > £i, L{t2, M) > £2, ■ " " , L{tm, M) > C) = det (1 + K^g) , (3.34) 



where the Fredholm determinant and Q{Tj,^) in the right hand side are defined in fl3.26p 
and fl3.27p respectively. The kernel is given by 

^3(n,6;T-2,6) 



where IC2 is the extended Airy kernel (13.28^ . 

This kernel has appeared in [SU |55] in the study of the height fluctuation property of 
the PNG model with external sources. When ri = ■ ■ ■ = = 0, it has been known that 
the Fredholm determinant of this kernel is described as the distribution of the larger of 
the largest eigenvalues in two independent GOEs, which is denoted as GOE^ [56l [57] . On 
the other hand, when ti = ■■■ = Tm = —00, the Fredholm determinant becomes the GUE 
Tracy- Widom distribution. Thus fl3.34p describes the transition of the largest eigenvalue 
distribution between GUE and GOE^. 

Next we consider the situation where the hopping rates qj {j = 1,2, ■■ ■ ,n) are more 
or less the same. In addition to the scaling fl3.2ip - fl3.24p . we set the hopping rate of the 
defect particles as 




(3.35) 



qi = q-q{l- q) 



D{u)mI 



1 5 



(3.36) 
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where D{u) is defined in (13.241) . Note that the parameters T]i{> 0){i = I,-- - ,n) char- 
acterize the inhomogeneity of the hopping rates and the case f]i = f]2 = ■ ■ ■ = Vn = ^ 
corresponds to the situation where the hopping rates degenerate completely. 
Under these settings, we obtain the following theorem. 

Theorem 2-3' 

lim Prob (L(ti, M) > i^, L{t2, M) > 4, ■ ■ ■ , ^(^m, M) > C) = det (1 + JC'^G) , (3.37) 

M^oo 

where the Fredholm determinant and Q{Tj,C,) in the right hand side are defined in fl3.26p 
and fl3.27p respectively and the kernel is given by 

"'If ( iw^\ ^ 1 

^3(n, ^i; T2, 6) =/C2(ri, ^i; T2, 6) + y] ^ / dwi exp i^iWi + ^)T\ — — 

^ 27r V 3 y % - n + twi 

1 f°° f iw'^\ 

X — / dw2 exp I i^2W2 + I [[{Vk -T2 + iw2), (3.38) 

where the contour runs from —00 to 00 passing the down side of the points «(% — Ti) {k = 
I,--- ,n). 

In the case m = 1, this Fredholm determinant has appeared in [l2l HSl [58]. Note that 
if we set rjj = for some j and 77^ = 00 for k 7^ j, we can realize the situation where 
the maximum of qi {i = 1,2, ■ ■ ■ ,n) is unique such as the former case. Indeed in the 
setting above, the kernel /C3 f l3.38p is reduced to /C3 fl3.35p . In that sense, the theorem 
is the generalization of Theorem 2-3. Later we will give only the proof of the theorem in 
Section 15. 2[ 

The Fredholm determinant fl3.37p also describes the limiting largest eigenvalue distri- 
bution in the following multi-matrix model. Let {-ffi}i=i,...,m be the N x N Hermitian 
matrices. Analogous to (13.300 . the joint density function of the model is defined as 

P {Hi, ti] ■ ■ ■ ; Hm, tm) 

= Z,.exp(-tri7o^)nexp (^^^^^^^l^^^^Z^^^^ exp{tTVHj. (3.39) 

Here Zy is the normalization constant and ti {i = 0, ■ ■ ■ , m) are parameters of the model 
such that t.j_i — ti gauges the strength of the connection between matrix Hi_i and Hi. 

V = diag(f 1, f2, ■ ■ ■ ,f„,0,0, ) represents the external source with rank n of the 

model. Under the scalings f l3.32p and 

V, = Vm (1-^], (3.40) 



the limiting joint distribution of the largest eigenvalues oi Hi {i = 1, ■ ■ ■ , m) is described 
as the Fredholm determinant (13.370 . 

lim Prob (/(^) < ai, ■ ■ ■ , /^""^ < ajj = det (1 + JC'^g) . (3.41) 



AT— >oo 
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In the case n = 1, fl3.4ip was shown in [55]. We can also show this equation for the general 
n case in a similar manner. 

3.2.5 Region 4 {u^ < u) 

This is the region where the effect of the defect particles has become dominant for the 
dynamics of the tagged particle. 

As in the region 3, we first consider the case 7^ ^2 7^ ■ ■ ■ 7^ ^n- In this case, one can 
expect that the dynamics is effectively the same as that of the slowest defect particle with 
stay rate q = max(gj). To see nontrivial correlations, times tj {j = 1, ■ ■ ■ ,m) should be 
macroscopically separated. We set 

ij = AGiuj)M - DGiuj)Mhj, (3.42) 

where uj = tj/M, Aciuj) is given in ( 13.13p . and 

Under this setting and introducing the parameter Tj such that 

e^^ = Dciuj), (3.44) 

we get the theorem as follows. 
Theorem 2-4 

lim Prob (L(ti, M) > h, L{t2, M) > £2, ■ " " , L{tm, M) > = det (1 + JCcG) , (3.45) 

M^oo 

where the Fredholm determinant and Q{Tj,C,) are defined in f l3.26p and (13.271) respectively 
and 

^ , forTi>T2. 
This kernel represents the propagation of a Brownian particle (a particle obeying the 
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Ornstein-Uhlenbeck process). Actually the two point distribution is calculated as 



lim Prob M) > £i, L{t2, M) > £2] 



+ 



1 

+ 2 



Si 



Sl 



d^2 
d^2 



K^n{T2,^2]n,C.l) /Cn(T2, 6; T2, 6) 

^n(^2,6; ^2,6) ^n(7'2,6;n,6) 



-(S2-'=^l~""2gl)^ 



rf6 



S2 g i_e2{n-^2) 



d£2 — , 
00 ^vr(l - e2(n-r2)) 



-f2 



TT 



(3.47) 



Next we consider the case gi ~ ^2 ~ ■ ■ ■ ~ ^n- The typical picture of the time evolution 
of the particles is depicted in Fig. 2(c). In this figure, it seems that four groups of which 
each top particle is the defect particle are formed and each group behaves like one particle 
in the TASEP. Thus in the region, one can expect that the system is effectively the same 
as n TASEP particles. 

In addition to the scaling fl3.42p - (13.441) . we also set 

qi = q-q{l-q)^. (3.48) 

M 2 



Note that the scaling exponent 1/2 in (13.421) is the same as that in the case of the central 
limit theorem. However the following theorem indicates that the limiting process of L{t, M) 
depends on the number n of the defect particles. Note that this is not the information 
which the scaling exponent can detect. 

Theorem 2-4' 

lim Prob (L(ti, M) > ii, L(t2, M) > ^2, ■ ■ ■ , ^(^m, M) > im) = det (l + /C^")^) . (3.49) 

M— >oo 



The Fredholm determinant andQ{Tj,^) in the right hand side are defined in (I3.26P and (13.271) 
respectively and the kernel is given by 



Ti,Xi;T2,X2 



) = ir*^")(ri,Xi;r2,X2) - 0(ri, Xi; r2, X2), 



(3.50) 



where 
K^''\n,^,;T2,^2) 



{2TTiy 



dwi / dw 



exp(w| -wf- 2w2^2 + 2wi^i) -A- e "^^^2 + e 



n 



e '^'^wi + ej 



exp 



(x2-en-^2gi)2 
I_e2{n--r2) 



for Ti < Ta, 

for Ti>T2. 



(3.51) 
(3.52) 
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In (I3.5ip . r represents the contour enclosing — e'^^e^ {i = I,-- - ,n) anticlockwise and 7 
represents the arbitrary path running from —ioo to ioo. 

The kernel above has been given in [55l [59] . As in Theorem 2-3', it is reduced to 
K-G (I3.46P if we set e = for some j and = 00 for A; 7^ j. Thus, later we will give only 
the proof of this theorem in Section 15.31 

The Fredholm determinant also describes the joint distribution function of the largest 
eigenvalue in the Hermitian multi- matrix model (13.391) . where we assume the rank of 
Hj (j = 1, ■ ■ ■ , m) is and V = diag(ei, €2, - ■ ■ , e„). The joint distribution of the largest 
eigenvalues Z^*-* of Hi can be described as the Fredholm determinant in the above theorem, 

Prob (/(^) < si, ■ ■ ■ , < sj) = det (1 + /C(")^) . (3.53) 

Note that in the equation above, we do not take the limit as the rank n goes to infinity. 

3.3 Scaling limit 2 (M fixed) 

Here we consider another scaling limit such that M is fixed and t goes to infinity. We set 

ti = e^^'T, (3.54) 
ii = il-q)ti-s,^2q{l-q)ti. (3.55) 




(3.56) 



and consider the asymptotic behavior as T goes to infinity. We get the following theorem. 
The proof of the theorem will be given in Section 15.41 

Theorem 3 

lim Prob (L(ti, M) > 4, L{t2, M) > £2, ■ ■ ■ , ^(^m, M) > C) = det (l + IC^^^^g) , (3.57) 
where the right hand side is the same as that of (I3.49p with n = M. 

3.4 Continuous limit 

In this subsection, we consider the continuous time version of the TASEP which is also 
usually discussed in the study of the ASEP. In the case without defect particles, the rule 
is defined as follows. Let t be the time variable which can take any positive real number. 
Between time t and t + dt, a particle can hop to its right neighboring site with probability 
dt if the site is empty. If the site is occupied, the particle stays at the same site with 
probability 1. 

We easily find that the continuous time version can be realized from the discrete time 
version defined in Section [2] by taking the limit, 

1 - g 0, t^oo, (3.58) 
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with t = (1 — q)t fixed. In addition, the scahng hmit in Section \32\ can also be studied by 
taking the t oo, M — > oo hmit with their ratio 

u = i/M, (3.59) 

fixed. Applying (13.581) and (13.591) to the Theorems in Section 13.21 we can rewrite the 
theorems for the continuous time TASEP. Here we show only the result which corresponds 
to Theorem 2-2. Analogous to (I3.2ip and (13.221) . we scale time tj and the position ij as 

ij = uM + C^''\u)Mhj, (3.60) 

ij = Ai'\uj)M - D^'\u)Mhi, (3.61) 

where 1 < tt and 
Then we have 

lim Prob (L(ti, M) > £i, L(t2, M) > 4, ■ ■ ■ , ^(L, M) > £j) = det (1 + /Cs^) . (3.65) 
Here the right hand side is defined in (I3.26p - (I3.28p . 

4 From the TASEP to the Schur process 

In this section, we discuss a relationship between the time evolution of a tagged particle 
in the TASEP with the step initial condition described in Section [2] and the stochastic 
growth of a Young diagram [521 E21 IMl ESj which turns out to be a special case of the Schur 
process jSHl Elj. In the last part of this section, the proof of Theorem 1 in Section [XT] is 
obtained. 

In the study of the TASEP, the technique using the enumeration of Young diagrams 
was used in [1^ to study the current fiuctuation. The discussion in [17] was based on 
the relation of the TASEP to the directed polymer problem on a matrix of which the 
elements obey the geometric distribution. In this section, on the other hand, we consider 
the correspondence between the directed polymer problem on a 01 matrix [28], [291 EQl [31] 
and the TASEP, which is a similar but different mapping. For the analysis of the tagged 
particle problem, the mapping is more natural and convenient. 



(v^-l) , 



M6 I a/S — 1 1 . 



(3.62) 
(3.63) 
(3.64) 
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4.1 TASEP and 01 matrix 



Let us consider an ensemble of iV x M 01 matrices {a{i, j)}i=i^... ^N,j=i,--- ,m, whose matrix 
elements are either or 1. All matrix elements are independent random variables and each 
element obeys the Bernoulli distribution, 

^^^^ ^ + 1 _ ^ |0, with probability 1 - g,, 
ll, with probability qj. 

Here 1 — qj is taken to be the same value as the hopping rate of the jth particle of the 
TASEP defined in Section |2l Each realization of the 01 matrix a{i,j) can be translated to 
a time evolution of the TASEP as follows. 

a{i, M + 1 — j) = 0(resp. 1): Between time i + j — 2 and i + j — 1, the jth 
particle tries to hop to the right neighboring site (resp. stays at the same site). 
When the right neighboring site is occupied at time i + j — 2, particle can not 
move due to the exclusion effect even if a{i, M + 1 — j) = 0. 

Since a{i, M + 1 — j) = could mean both hopping or stay depending on whether the 
target site is occupied or empty, the mapping from a 01 matrix to the time evolution of 
the TASEP is not one-to-one. In Fig. 4, we depict an example of the time evolution in the 
TASEP and the corresponding two 01 matrices. Note that in the step initial condition, the 
jth particle cannot move until time j — 1 and that the time evolution of the jth particle at 
time step between i and i + 1 does not depend on that of the j — 1th particle at the time 
step. Thus one finds that the 01 matrix defined above has all informations needed for the 
time evolution of the Mth particle from time through M + N — 1. 

Let us define a left-down path tt{N,M) of a N x M 01 matrix {a{i, j)}i=i,... ,N,j=i,--- ,m 

by 

7r{N, M) = {{{ik, Jk)}i<k\a{ik,jk) = lA < < ^2 < ■ ■ ■ < N, M > ji > j2 > ■ ■ ■ > 1} . 

(4.2) 

Note that the row indices ik {k = 1, ■ ■ ■) are strictly increasing while the column indices 
jk {k = 1, ■ ■ ■) are weakly decreasing. In vr, the position of {ik+i,jk+i) is always located on 
the down side or left-down side of the position of {ik,jk) in the 01 matrix. We also define 
the quantity G{N, M) by 

G{N,M)= max WiN,M)l (4.3) 

tt{N,M) 

where |7r| denotes the number of elements in the path vr. If we regard the 01 matrix as 
the random media in a plane, 7r(A^, M) as a spatial configuration of a polymer chain, and 
|7r(A^, M)| as its energy, we can interpret (14.31) as a statistical mechanical problem of a 
directed polymer. This kind of problem is studied in [28l [29l [301 EI]- 

Now we find that this maximum length of the left-down path is directly related to the 
position of a tagged particle as 

G{N, M) = d{N, M) (4.4) 
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where 

d{N,M) := N - L{t = N + M -1,M) (4.5) 

represents the number of times the Mth particle stays at the same site from time through 
N + M — 1. The equation (I4.4p enables us to study the dynamics of the TASEP in terms 
of the 01 matrices and plays a fundamental role in our subsequent analysis. 

We can show (14. 4p in the following way. For the case M = 1 (the case where we tag the 
first particle), we can easily check this equation, since G{N,1) defined in (14. 3p is simply 
the total number of the value 1 in the sequence {a{k, l)}k=i,...,N and it is clear that this 
quantity represents the number of times the particle stays from t = through N. 

For the case M > 2, we can prove (14.41) by mathematical induction about the row A^. 
Here, we mainly consider only the case M = 2. For the case M > 3, we can check (14.41) in 
a similar fashion to this case. 

When = 1, we can check ( 14.41) easily for all four cases (a(l, 1) = 0, 1 and a(l, 2) = 
0,1). 

We assume that (10|) holds for A^ = Na, and M = 1 and 2, i.e., 

d{NaA)=G{NaA), d{Na,2)=G{Na,2). (4.6) 

Under this assumption, we will show that d{Na + l, 2) = G{Na + l, 2) as follows. From (14.31) . 
one easily finds 

G{Na,2)>G{K,l), GiNa,t)>G{N,,t), (4.7) 

where 1 < A^^ < A^^, and i = 1 and 2. Considering these properties, we classify the problem 
into two cases, 

(i) 

G{Na,2) = G{NaA), and a(Ar, + 1, 2) = 1, (4.8) 

(ii) 

^(A^,, 2) > GiNa, 1) and a{N + 1, 2) = or 1, 
G{Na,2) = G{NaA) and a(A^ + l,2) = 0. ^ ' ' 

The case (i) corresponds to the situation where the two particles occupy the neighboring 
sites at t = A^a + 1- Thus the second particle (the tagged particle) must stay at time step 
between A"a + 1 and A'a + 2 . Hence one finds 

diNa + l,2) = diK,2) + l, (4.10) 

regardless of the value of a{Na + 1, 1). On the other hand, due to (14.31) . (14.71) and (14.81) . 
we have 

GiNa + 1, 2) = GiNa, 1) + a{Na + 1, 2) = GiNa, 2) + 1. (4.11) 
Thus from (14.61) . (I4.10p and (14.111) . we have for the case (i), 

diNa + l,2)=GiNa + l,2). (4.12) 
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The case (ii) corresponds to the situation where at t = Na + 1, the distance between 
the first and the second particles is at least one site. Note that in this case we do not need 
to consider the exclusion effect of the first particle and thus the situation is essentially the 
same as the case M = 1. We have 

d{Na + 1,2) = d{Na, 2) + aiNa + 1,1). (4.13) 

On the other hand, from (14. 7p and (14.91) we can also check that 

G{Na + 1,2) = G{Na, 2) + a{Na + 1,1). (4.14) 

Thus we find that (I4.12p holds also for the case (ii) and therefore (14. 4p holds for M = 2. 

In the case M = k {k > 3), we can also classify the situation into two cases as in (i) 
and (ii) for M = 2, i.e., the case (i)' where all k particles are packed and (ii)' where the 
particles form several groups which are separated by some successive empty sites. For the 
case (i)', we can check (14.41) in a similar manner to the case M = 2. For the case (ii)', we 
can easily show that the case is essentially the same as the case (i)' for M = m which is 
smaller than k. 

4.2 01 matrix and the dual Robinson-Schensted-Knuth corre- 
spondence 

The quantity we would like to discuss is the multi-time distribution of L{t, M), 

ProbTASEp(^(ti, M) > £i, L{h, M) > £2, ■ ■ ■ , Htm, M) > £„), (4.15) 

where ProbTASEP represents the probability measure of the TASEP defined in Section [2l 
Due to the result in Section l4?]T especially (14. 4p and (14. Sp ). one finds 

ProbTASEp(^(iVi + M - 1, M) > £1, 

L(iV2 + M - 1, M) > £2, ■ ■ ■ , L{Nm + M - 1, M) > £^) 
= Proboi(G(iVi,M) <iVi-£i, G(iV2,M) <iV2-£2,--- , G(iV^, M) < iV^ - £„), 

(4.16) 

where Proboi represents the probability measure of 01 table defined in (14. ip and G(iVj, M) 
is defined in (14.40 . 

In this section, we further restate our problem in terms of the combinatorics of Young 
tableaux. For the definitions and basic properties of Young tableaux and related subjects, 
we refer the readers to [6OJ. There is a bijective mapping between an x M 01 matrix and 
a pair (P, Q) where P* (the transpose of P) and Q are the semistandard Young tableaux 
(SSYT) with the condition that the shape of P is the same as that of Q. The mapping 
is called the dual Robinson-Schensted-Knuth (RSK) algorithm. In order to get {P,Q), we 
first construct a two-line array (or generalized permutation), 

/ 1 ... 1 2 ■■■ 2 ■■■ ■■■ N 

( ,(1) . . . ,(2) . . . ,(2) . . . AN) . . . AN) 

\Jl Jmi J I Jm.2 J I Jmjv 



(4.17) 
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where 1 < j[ < j2 ■ ■ ■ < M, from the 01 matrix by hsting (*.)'s for which a{i,j) = 1. 
Next, we construct P (resp. Q) by arranging on a plane the figures in the second row 
(resp. the first row). Note that each positive integer constructing P(resp. Q) is less than 
M(resp. N). For details of the algorithm, see [60j. The generalized permutation and 
(P, Q) corresponding to the 01 table in Fig. 5(b) are given in Fig. 6. 

Now we express G{N, M) fl4.3l) as the quantity related to a pair (P, Q). Let us denote 
the shape of P (or Q) by the Young diagram, \{N, M) = {\i{N, M), A2(iV, M), ■ ■ ■ ) where 
Aj > represents the length of ith row of the tableaux P (or Q). In the example of 
Fig. 6(b), we see A = (4, 3, 2, 2, 2). Then we have 

G{N,M) = X[{N,M). (4.18) 

Here A' = (A'^, Ag, . . . ) means the transpose of the Young diagram A. Thus X'i{N, M) is 
equal to the length of the first column of A(A^, M). 

Eq. f l4.18p may be shown by examining the dual RSK algorithm directly but is also 
understood as follows. For the sequence of the numbers {ii,i2, ■ ■ ■ ) where i G (0, 1, 2, ■ ■ ■ ), 
we define a nondecreasing (resp. nonincreasing) subsequence (ji,j2,''') such that ji < 
32 < ■■■ (resp. ji > j2 > •■■)• From the construction it is not difficult to see that 
G{N, M) is equivalent to the length of the longest nonincreasing subsequence of the second 
row ■ ■ ■ ,jml, ■ ■ ■ ■ ■ ■ , jS) ill a generalized permutation fl4.17p . Thus, G{N, M) 

is also regarded as the length of the longest nondecreasing subsequence of the opposite 
sequence {jin2, ■ ■ ■ ■ ■ ■ ,jml, ■ ■ ■ jJi'^) of the second row in fl4.17p . We can apply the 

(normal) RSK algorithm JEU\ to this reverse sequence to obtain a SSYT. Let us denote 
the SSYT by R. Then ( 14.181) is a consequence of the facts that the length of the longest 
nondecreasing sequence in the reverse sequence is equal to the length of the first row of 
R (this is known as a property of the normal RSK algorithm [171 120]) and the sjTumetry 
property, 

P* = P, (4.19) 

whose proof can be found in Appendix A of [61]. Notice that in the example of Fig. 5(a), 
L(9,4) = 1 thus rf(6,4) = 5 while in Fig. 6(b), A; = 5. 
Thus from fl4.18p . we get 

Proboi(G(iVi, M)<Ni-iu G{N2, M) < N2 - £2, ■ ■ ■ , G{N^, M)<N^- i^) 

= FiohxiX[iN,,M)<N,-i,, \\{N2,M)<N2-l2r-- ,K{Nm,M)<Nm-U. 

(4.20) 

Here ProbA in the right hand side represents the probability measure of the set of Young 
diagrams obtained by the dual RSK algorithm from the 01 matrix defined in (14. ip and 
A(iVj, M) means the Young diagram obtained by applying the dual RSK algorithm to the 
A^i X M 01 submatrix {a{i, j)}i<i<Ni,i<j<M of the A^ x M 01 matrix (14. ip . 

From (14.161) and (I4.20p . we managed to interpret the time evolution of a tagged particle 
in the TASEP as the time evolution of A'^, the length of the first column, regarding Aj as 
a time parameter. 
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4.3 Growth of Young diagrams and Schur process 



In order to investigate the right hand side of (]4.20p . we consider the following probability, 
ProbA(A(l, M) = AW, A(2, M) = A^^), ■ ■ ■ , X{N, M) = A^^^ = A) (4.21) 

for a given set of Young diagrams, =i,...,N- This equation describes the growth of 

a Young diagram. Eq. fl4.2ip represents the joint distribution of them. Fig. 7 illustrates 
the growth of Young diagram which corresponds to Fig. 5(b). Notice that from (14. 4p 
and (14.181) . the growth of X'i{i,M) describes the time evolution of Mth particle in the 
TASEP. Let us denote the Q^*"'-* as the part of the Q tableau obtained by picking up only 
the symbols 1, 2, ■ ■ ■ ,i. From the rule of the dual RSK algorithm [60], we find that Q*-'^-* 
is equivalent to the Q tableau obtained by applying the dual RSK algorithm to the the 
submatrix {a{i, j)i=i^,,,^k,j=i,...,Ai}- Thus one has 

sh(QW) = A(i,M), (4.22) 

where sh((5''*^) denotes the shape of Q^^\ Hence the Q tableau records the growth of 
Young diagram, that is, it one-to-one corresponds to the set of the Young diagrams 
{A(z, M)}j=i ... AT. Compare the tableau Q in Fig. 6(b) with {A(i, M)}j=i_... g in Fig. 7. 
We can recognize that sh.{Q^^^) constructed from the tableau Q in Fig. 6(b) is equivalent 
to X{i,M) in Fig. 7. Hence we have 

ProbA(A(l, M) = A(^), A(2, M) = A^^), ■ ■ ■ , X{N, M) = A^^)) 

= Prob|p,Q}(sh(Q('^)) = A«, A; = 1, . . . , iV). (4.23) 

Here Probjp^g} represents the probability measure on the set of the pairs {P, Q} obtained 
by the dual RSK algorithm from the x M 01 matrices defined in (14. ip . Hence it is 
found that the growth process of the Young diagram is characterized by the Q tableau 
in the dual RSK correspondence. This type of growth process has been also discussed 
in [321 [331 [31 [35]. 

Remembering P is constructed from the second row of a generalized permutation, we 
find that the symbol j in P represents the column index of the element of the 01 matrix 
where a{i,j)=l. Thus the total number of js in P is equal to the total number of the 
figure 1 in the jth column of the 01 matrix. Hence, for a given P', one finds 

Prob|p,Q}(P = P', sh(Q('=)) = X^'\k = 1, ■ ■ ■ , iV)) 

(<lTW? ■ ■■<ir''' X (1 - gM)^-«^(^')(l - gM-i)^-«^(^') ■ ■ ■ (1 - qr-'''(^'\ 
= < if A*-"*^^ and A'-^^^-'/A^*^ (i = 1, ■ ■ ■ , — 1) have no two squares in the same column, 
[O, otherwise. 

(4.24) 

where tl^(-P') means the total number of the symbol i in P' . Applying the above equation 
to the combinatorial definition of the Schur function, 

s 
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where the summation is taken to all SSYT S with shape A//i, we have 

Prob{p,Q}(sh(Q('=)) = AC^), k=l,...,N) 

= Prob|p,Q|(P = P\ sh(Q«) = AW, A; = 1, . . . , iV) 
p' 

= (1 - qi)^ ■■■{!- qM)'^Sy(N){pM, ■■■ 

s;,(i)(l,0, ■ ■ ■)sa(2)/a(i)(1,0, ■ ■ ■) X ■■■ X SA(iv)/A(iv-i)(l,0, ■■■), (4.26) 

where Pi = qi/{l — qi). The factor Sa(i)(1,0, ■■■)■■■ ensures that the probability is zero 
unless AW and A^^+^V^^*^ = 1, ■ ■ ■ , - 1) have no two squares in the same column. 
Thus combining (14.231) with (14.261) . we finally find 

ProbA(A(l, M) = AW, A(2, M) = A^^), ■ ■ ■ , X{N, M) = A^^)) 

= (l - qi)^ ■ ■ ■ {I - qu)'^ Sy(N){pMr ■ ■ ,Pi)x 

Sa{i)(1,0, ■■■)sA{2)/Aa)(l,0,---) X ■■• X Saw/;^(^-i)(1,0,- ■■)• (4.27) 

This measure is the special case of the joint distribution function in the Schur pro- 
cess [Ml ST]. Since the Schur function is expressed as a determinant, the function (I4.27P 
is described as a product of determinants. Furthermore, it is known that the process has 
a remarkable mathematical structure that the correlation function of {X[''^}i=i,2,--- , j=i,---N 
can be expressed as a determinant and its scaling limit can be discussed exactly. The Schur 
process also appears in other fields of physics and mathematics such as random growth 
process [38], melting problem of a three-dimensional crystal [62], random tiling model [63] . 
and so on. 



4.4 Proof of Theorem 1 

From (|4A6|) . (14:201) and KTl} . we find 

ProbTASEp(^(ti, M) > ^1, L(t2, M) > £2, ■ ■ ■ , L{t^, M) > £„) 

= Probsp(A;(iVi, M) < iVi - A^iVa, M) < - £2, ■ ■ ■ , KiN^, M) < - C). 

(4.28) 

Here Probsp represents the probabilistic measure of the Schur process (14.261) . 

In general, the Schur process is defined as follows. For the set of the Young diagrams 

A(^) D /x(^) C A(2) D ■ ■ ■ D C A(^), (4.29) 

the following weight is assigned 

Sxm (Po )sa(i) /^(i) ipl)sx(2)/^,w (pt) ■ ■ ■ Sx(n)/^(n-i) (p^_i)sA(iv) (p^). (4.30) 

Here Sp/^{p) is the Schur function with the specialization of algebra p. In our case, we 
notice only two cases 

.A/» = |'^/^^"^'---'""^' f-^ = /^^^) = n|iT3fe> ^4.31^ 
[sy/^Ao-u ■ ■ ■ ,ajv), for p = p{z) = Uj=ii^ + S^)' 

22 



where sx/^{ai, ■ ■ ■ , oat) is defined in (14.251) . Substituting to (14.291) . 



pr{z) = l + z, {l<^<N-l), 




(4.32) 



where Pi = qi/ (1 — g^) and applying Theorem 1 in [36] or Theorem 2.2 in [37] to this case, 
we finally obtain the Fredholm determinant representation( (I3.3P - (13. 71) ) of (14.281) . 

5 Asymptotics 

In this section, we discuss the scaling limit of the multi-time distribution function (12.11) by 
applying the saddle point method to the kernel (I3.5l) - (l3.7p in Theorem 1. 

We consider two types of scaling limit. The first case is explained in Section [221 We 
focus on the situation where both time t and the label of the tagged particle M tend to 
infinity when there are n defect particles with the stay rates „ in front of the 

tagged particle. The results on the limiting distributions are summarized as Theorem 2. 
Their proofs are given in Section 5.1, 5.2 and 5.3. 

In the second case, we take t to be infinite with M fixed. The limiting distribution is 
described in Theorem 3 in Section 13.31 whose proof is given in Section 5.4. 

5.1 Region 1 (proof of Theorem 2-1) 

In this section, we consider the scaling limit in the region 1 in Fig. 3. We scale time ti 
as (I3.14p . In order to calculate the limiting distribution, we rewrite the kernel (I3.5l) - (l3.7p 
with the condition (13. 9p as 



1 f dz, r dZ2 Z, (1 + l/z2Y^-''+' ( 1-PZ2 \ A 1 - P^Z2 

{2m)^Jn, Z, Jc. Z2 ^1 - Z2 (1 + lM)*i-^^+l VI J \\^-V^ZxZr 

(5.1) 



where p = 5'/(l — (?) and "pi = qi/{l ~ Qi)- Cr denotes the contour with radius R enclosing 
the origin anticlockwise and R[ {i = 1,2) satisfy the condition, 1 < i?2 < -R'l < 1/pi for 
ti > t2 and 1 < R[ < R2 < 1 / pi for ti < ^2- Using (15.11) and the relation 




(5.2) 
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one finds 



A-^f ^-i ^ ^ _ )T.m=o^i{xi+^M)^2{x2 + mM). ti>t2, 
- Lm=o^u^i - - l,ti)^2(a;2 - m - l,t2), h < t2. 



Here the function \l/i(m, ti, xi) and \E'2(^5 ^2, 2:2) are given by 

1 ( . 1 f dz 1 -A- 1—^2; 1 

Mx2,t2) = ^f ^{l + l/wy^-'^^\l-pwrfl\^w-\ (5.5) 



wliere tlie radius Ri of tlie contour C/?^ in (15.41) is taken such that 1 < i?i < 1/pi. In the 
following discussion, we evaluate the asymptotics of the kernel (15. ip by applying the saddle 
point method to the functions \E'j(x,t) (i = 1,2). 

First we consider the asymptotics of \E'i(xi,ti). Substituting (I3.14p with (15.40 . we set 

M^ut.) = ^f ^e-««') L^fl ^-^:- J (5.6) 

where 

/(z) =plog(l + 2;)+log(l-pz). (5.7) 
From the equation f'{zc) = 0, we find the critical point Zc is given by 

z, = 0. (5.8) 
Scaling the variable z around the critical point as 

z' z' 

z = z^^ r = r, (5-9) 

D1M2 D1M2 

where Di = y/q/{l ~ q) = a/p(1 + p), we evaluate the asymptotics of the term e~^^'^^'^ 
in (15. 6p by the saddle point method. 



-Mf{z) ^ y'' ' 2d(m J 



,'2 



e-'^'f^'>^e V" ^i^i^' ;=e-, (5.10) 

as M — s> 00. Next we consider the asymptotics of other terms in (15. 6p . From (15.90 . one 
obtains 

1 " 1 

5 



nf^~i. (5.11) 

-L-L 1 - ViZ 



Transforming Xi into x'j {i = 1,2), 

= ti -M + 1 (5.12) 
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and considering (15.91) . one gets 

1 T'^'l 

(5.13) 



^xi-ti-\-AI-l 



1 



Note that by the transformation fl5.12p . the function g{ti,x) (13.40 in the Fredholm de- 
terminant (13. 3p is changed to giji^x') = —Xi~ooA)i^')- From (I5.9I) - (I5.13I) . we obtain the 
asymptotic form of \l/i(xi,ti), 

^^^Dr]M^ r"^^'e4-.^-,-.;-. (5.14) 



ico+e 



where e > 0. 

Next we consider the asymptotics of '^2{x2,t2) (15. 5p . Similar to \l/i(xi,ti), we set 

vi/2(a;2,t2) = — /" ^e^^("')(l + ^)i+^*^i^^fT^^^^x;-^-*+^-^ (5.15) 
2Tci J(j^ w 1 

Here f{w) is defined in (15. 7p and Ci denotes the contour enclosing the origin anticlockwise 
with radius 1. Applying the saddle point method to this equation in the same way as 
^i(a;i,ti), we get 



^y^^_,-^..». (5.16) 

Thus from ([53D, dSH, and (IHTTBD . we find 
lim K{ti,Xi;t2,X2) 

-(Z}lV^)-2-<E"=oV^lK+"^+l'^l)^2(4+^+l,r2), n<T2, 

where 



(5.17) 



V'i(x;,ri) = — / dze'^-^^'z<-\ (5.18) 

Mx'2,r2) = ^ [ 4^e-^--. (5.19) 
Introducing the function Dn{z) by 

D,_i(r) = v^e"^V^i(x, r), (5.20) 
one finds that the function satisfies Weber's equation |48j . 

'^'^"^''^ ^n + l- -z''^ D^{z) = 0. (5.21) 



dz'^ V 2 4 
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The initial condition was given in [30] as 

D„(0) = sin [^{n + 1)) r j , (5.22) 

D;(0) = -Z^„+i(0), (5.23) 

where T{x) represents the Gamma function. The function satisfying (I5.21l) - fl5.23l) is known 
as the parabohc cyhnder function [48J. In addition, one easily finds il'2{x, t) can be repre- 
sented as (13.181) from the integral representation of the Hermite polynomial Hn{x) [48j, 

'^»(^) - / iSr^^"-. (5.24) 



Hence, noticing the prefactor (Div^M)^2 ^'i does not contribute to the determinant 
in 03.31) . we get the limiting kernel (I3.16l) - fl3.18l) . 

5.2 Regions 2 and 3 (proof of Theorems 2-2 and 2-3') 

In this section we consider the scaling limit of the regions 2 and 3 in Fig. 3. We fix the 
scaled time u (13. Sp such that 1/(1 — g) < u < Uc for the region 2 and u = Uc for the region 
3, where Uc = {(f — 2qq + q)/{q — qY- We scale the time tj, the position of the tagged 
particle at tj, and the stay rates of the defect particles (z = 1, ■ ■ ■ , n) as (I3.2ip . (I3.22p 
and (I3.36P respectively. 

First we analyze the asymptotics of K{ti, xi, t2, X2) in (13.60 by the saddle point method. 
Using the variable fi{u) which will be fixed later, we set 

A(ti,xi;t2,X2) = ,^ / / exp (M(/„2(2;2) - /„i(2;i))) 



{l-Kif Jc^^ Zi Jc^^ Z2 ' ' 2:1-^2 



X 



1 - J92i V -p-r 1 - PiZ2 I + I/Z2 



n^^^^^^^. (5.25) 



\l-pz2J fj^ 1 - PiZl 1 + 1/-Z1 ' 

Here ui = U/M, p = q/ {1 — q), and pi = qi/{l — (ji)- The function fu{z) is defined as 

Uz) = iu-l) log(l + z)+ log(l - pz) + (/i(M) -u + 1) log(^). (5.26) 

We fix the value of /x(m) in such a way that two saddle points of fu{z) merge to one 
point. We have 



. / X p(m - 2) + 2x/p(u - 1) 

K^) = u-l- A2iu) = ^ -r-^ 5.27 

l+p 



and the double critical point Zc{u) of fu{z)is obtained as 



= -y^'-VP . (5.28, 

/p'^{u - I) + y/p 
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Note that 

fLi^ciu)) = = (5.29) 

is satisfied with fi{u) fixed as f l5.27p . 

Thanks to the property that Zc{u) < 1/pi when 1/(1 — g) < u < Uc which is the 
condition of regions 2 and 3, it is found that we can deform the contours Cr. [i = 1,2) 
around the double saddle point ( I5.28p . Thus we scale zi and Z2 as 

zi = Zc{ui) 1 — — - ] , Z2 = Zc{u2) 1 H — — r , (5.30) 

\ D{u)M^ J \ D{u)Mz ) 

where D{u) is defined in (13 .24^ . Noting (15.291) and the relation 

- rW|MM . o(uf, (5.31) 

we have 

/// 

, , , , exp{MU,{zM) + Mf-^^^{Z2 - zMf) 

exp(M/„,(z,(«i)) + M ^"^^f"^» (^i - z,(ni))3) 

~ exp (m(/„,(z,(«2)) - fuAzciu,))) + Ml^'^ . (5.32) 

Hence we had the asymptotic form of the term exp {M{fu2{z2) — /^^(-Zi))) in (I5.25p . Next 
we consider that of other terms in (I5.25p . Combining (15.300 with the relation 

Zc{u) D{u)M3 

we find 

Zi ~ Zc(u) 1 H r , -^2 ~ zJu) 1 H r • 5.34 

^ ^ \ Diu)Ml J ^ ^\ D{u)M\ ) ^ ' 

Furthermore, associated with (13.221) . we scale Xi {i = 1,2) as 

Xi = {u- 1)M - A2iui)M + D{u)M^ii. (5.35) 

Here Ai^u) is defined in (IXT^ . Using (IHT^ and ^M), we get 

zi D{u)Mh 



Zi - Z2 n- T2- i{Wi + W2) ' 



(5.36) 



l-pz,Yl + l/z2 _ ^ ^5 3^^ 



J92;2y 1 + 1/zi 

fri^./l. region 2 = tM?) , 

i|l-p.^i region 3 (« = i^J, 

~ z,(n)^(")^^«-«^) exp(6r2 - + + ^^26)- (5-39) 



x-i_—fi{u-i)M 
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Note that in fl5.38p . the asymptotic form is different between regions 2 and 3. This leads 
to the difference of the hmiting distribution between Theorems 2-2 and 2-3'. 
Thus we obtain the asymptotic form of K{ti,xi; t2, X2), 



K{ti,xi;t2,X2) 

dAe-^(---)Ai(6 + A)Ai(6 + A), 

JJ[u)M 3 

r /o°° dAe^^(-^--^)Ai(a + A)Ai(6 + A) 



region 2, 



D(u)M7 



+ 



D{u)M^ 



T E"=i h n [liiwi + ^ riLi 



1 



(5.40) 



^ ^ h I-oc '^^2 exp (^iC2W2 + ^) IliJiiVk -T2 + iw2), regiou 3. 



Here 



u = exp {M{U,{z,{u2)) - fuA^c{u^))) + - ^1^1) 
In (15.40 p . we used the integral representation of the Airy function 



Ai(x) 



2ti 



(5.41) 



(5.42) 



and the relation 



1 



T1-T2- i{wi + W2) 
1 



N 

n 



rjj -T2- IW2 

\ Vj - n + iwi 

r]i-T2- iw2 



^71 - n + iwi {rji - Ti + iwi) {ri2 - n + iwi) 

{rji -T2- iw2){r]2 -T2- iw2) ■ ■ ■ ivN-i -T2- iW2) 



+ 



(^71 - ri + iwi){7]2 -Ti+ iwi) ■ ■ ■ {r]N - n + iwi) 



(5.43) 



Next we consider the asymptotics of (j)ti,t2{^i^ ^2) (13.70 . Using fu{z) (15.261) and /i(u) (I5.27p . 
the function is rewritten as for ti < t2, 



2m 



dz 



t2(Xl,X2) = — — / —Q'^'^Uu2{^)-fu^{z))^X2~lJ.{u2)M-Xi+lJL{ui)M^ 



(5.44) 



Ci 



Let the variable z scale as 



Z = Zc{u) 1 + 



D{u)M 



1 I 1 

3 



where Zc{u) defined in (I5.28p is the double saddle point of ([51 
for z = 1, 2 



z ~ zJui) 1 + 



D{u)mI 



(5.45) 

Due to (15.331) . we have 
(5.46) 
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From the above two equations, (15.311) and fl5.35p . we get 



gM(/„2(2c(«2))-/ui(2c(Mi)))g-(r2-n)fT2-i(r|-r2)<j+J^^ ^5_47^ 



^^^^)Z.i«)M3(5,-5,,g..(5,-5,,_ (5.48) 
Thus we obtain for ti < t2, 

4>ti,t2i^i^^2) 

33/* 00 

_g€lTi-6T2 + ^S^ / ^^Q-{r2-riW-i{Tl-T2)a+i{i2-ilh 

u)MI J -00 



D{u)Mi 



V 



exp( -fe-eO;_fe + 6)(r.-r,)^(r.-n. ^^^^^^ 



D(m)M^ ^47r(r2 - n) V 4(r2 - n) 2 12 

where is given in (I5.4ip . From (I5.42p we finally find 

0ii,i2(a;i,^2) ~ — — ^ / rfAe-^(^^-^^)Ai(ei + A)Ai(6 + A). (5.50) 
D{u)M3 J^oo 

Hence from (15.430 and (I5.50p and noting that the term u (I5.4ip does not affect the deter- 
minant, we finally obtain the desired kernels for the regions 2 and 3. 

5.3 Region 4 (proof of Theorem 2-4') 

In this section, we discuss the asymptotics of the region 4 in Fig. 3. In this region we 
take time tj in such a way that they are macroscopically separated and the scaled time 
Uj = tj/M is taken as Uj > Uc- The position ij of the tagged particle at tj and the stay rates 
gj (z = 1, • ■ ■ , n) of the defect particles are also scaled as (I3.42p and (I3.48P respectively. 

Deforming the contour Cr-^ in the kernel K{ti,xi;t2,X2) (13. 6p . we divide it into two 
parts, 

k{ti, xi; t2, X2) = Kiih, xi; ts, ^2) - ^2(^1, t2, X2). (5.51) 

Here for i = 1,2 

ki{ti,Xi]t2,X2) 



dz, f dZ2 Z, (1 + 1/^2)*^-""+^ / 1 - p^2 \ fr 1 - P^Z2 ^2 



(27rz)2 J 0^2 ^2 ^1 -^2(1 + l/z,y^-''+^ \l - pzj H 1 _ p^^^ 



v^2 



"1 

(5.52) 

where p = q/{l — q) and pi = qi/{l — qi). In the equation above, the contour F^^"* encloses 
z = 1/pi anticlockwise while the contour F^^"* is chosen in such a way that it encloses 
z = —1,0,1/pi anticlockwise. 
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First, we discuss the asymptotic form of Ki{ti,xi;t2,X2)- This is rewritten as 
{2my JrW zi Jc^,^ zi - Z2 

1 - P^l y TT 1 - PiZ2 1 + 1/Z2 z^^'^^^^""^^ 

where A(-u) is a variable which will be fixed later and 

gu{z) = {u-l) log(l + 2) + log(l - pz) + (A(m) -u + l) log 2;. (5.54) 

We choose the value of A('u) in such a way that the saddle point of gu{z) is 1/p where 
p = q/{l — q) ■ From the condition, 

9u Q) = 0, (5.55) 

we have 

A{u) = {u-l)- Ag{u) = + (5.56) 

I + p p — p 

Scaling Zi {i = 1, 2) around the saddle point (I5.55P as 

z^='-(l- . V (5.57) 

and noting p^D'q{u) = 2g"(l/p), we can obtain the asymptotic form of e^^*^^"2(^2)-3ui(zi)) 
in fl5.53p by the saddle point method, 

^M{gu2{z2)-gu^{zi)) ^ gA^(5-2(i)-9«i(i))e^i-^f. (5.58) 
Associated with (13.421) . we set 

Xi = {ui - 1)M - Ag{u,)M + DG{ui)M'^i,, (5.59) 
for z = 1, 2. From this equation, (15.571) and (13.481) . we find 

n l -PiW ^ -A- e"^2^2 + Q zi Mi 



.^^ - PiZ e '^l^^;l + e^' - 2:2 -2(e '^^Wi-e '^^■012) 

l-pZ^Yl + I/Z2 4-MA(..) (e^25.-ena)M^ 



1~PZ2J 1 + ' 

where we used the parameter Xj defined in (13.441) . 



^_2^252+it«ia^ (5.60) 
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Thus we eventually get 



ii±Z^ P / dwi / dw2- n (5.61) 



Here the contour F encloses — e^^e^ (i = 1, ■ ■ ■ , n) anticlockwise and 7 is an arbitrary path 
running from —ioo to too. 

Next we consider the asymptotics of 0j^^j2(xi, X2) (13. 7p . By use of gu{z) (I5.54p and 
A(u) (I5.56p . we set for ti < t2, 



27rz ./ci ^ 



Scaling z as 



1 / 2u; , , , 

z = -\ l+i , 5.63 

v/(/^^(«,)-D^(«o)Mr 



From this equation and (15.590 . we get 

gM(9„2(^)-g„,(^)) _ gA^(9..2(i)-9-i(|))e-«''^ (5.64) 



,X2-xi-M(A(«2)-A(ni)) _ / i 1 ^^r, o,'.,. ^ 



Thus we obtain for ti < t2, 

^,.,M,X2)r^^^.^== dwe V^^^ 



(^1 _|_ py^^^ pXl-X2 g I_e2(ri-r2) 

e^2M^ A/7r(l - e2(n-r2))' 



(5.66) 



At last, we discuss that K2(ti,xi;t2,X2) in (I5.52p does not contribute the asymptotic 
form of the kernel (13.51) . The derivation can be done in the similar fashion to the one in 
Theorem 3.1. in [M]- Here we give only its outline. 

First we consider the case of equal time, ti = t2 = t = uM. Scaling Zi {i = 1, 2) 
as (I5.30p . we have 



K2{t,xi;t,X2) 



r\-/ 



D{u)M3 

fAA(u) 2 Dciu) 1 \ . fAA(u) 2 DgM 1 
V D{u) D{u) J V D{u) D{u) 

(5.67) 
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where AA{u) = A2{u) — Ac{u). Noticing that AA{u) > for the region 4 {uc < u) and 
the asymptotic form of the Airy function 

Ai(x) ~ ^ ^ x~^ exp ( — ^a^M , (5.68) 



2v^ "V 3 
as X — » oo, we find 

K2(t,xi;t,X2) -e-'^^^'). (5.69) 

This indicates that under the scahng fl3.42p . the kernel K2{t,xi;t,X2) vanishes as M goes 
to infinity. 

Next we consider the case for arbitrary ti {i = 1,2). We scale Zi around the double 
saddle point Zc{ui) of fmiz) and Z2 around the saddle point 1/p of (^ujl-^i) as 

z. = z.{u,)\l--^],z. = Ul (5.70) 

Under these scalings, we find 

K2{t,x,-,t,X2) - ,\ ^^^^^^ X e(«i,ei), (5.71) 



where 



e(«i,6) 



1 ( l-p/p "-Q_l-Pi/P ZciUi) ""^p 

1=1 ^ 



Zc{ui) -1/p \1 -pzc{ui)j fj^ l-piZc{ui) D{ui)Dg{ui)M-6P 

Considering the case ti = t2 in (15.711) and the former result (15.691) . one easily finds 

e(M,0 ~e-^(*'). (5.73) 

Hence we finally get 

i^2(ti,xi;t2,X2) ~ {I+pY'-'^p^^-^^ (e-^(*^)) . (5.74) 

From (15.611) . (15.661) and (15.741) and noting the factor (1 -\- py2~tip^i-^2 fjoes not affect 
the determinant, we obtain the desired expression for the limiting kernel. 
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5.4 Fixed M case (proof of Theorem 3) 

In this section we discuss the scahng hmit explained in Section 13.31 where time t goes to 
infinity with M fixed. In this limit, we scale the time tj, the particle position £j, and the 
stay rates (i = 1, ■ ■ ■ , M) as (13.541) . fl3.55p . and fl3.56p respectively. 

First we consider the scaling form of K{ti,xi\t2.,X2) (13. 6p . Changing the variables 
Zi {i = 1, 2) to z'^ = l/zi, we have 

K(t x-t x) = ^ f ^-^ f ^4 4 (i + 4r--^-^+^ fT 4-i>.4"^^ 
^ " " " i^my 4 z[ Jc,^ 4 4 - (1 + z[y^-M+^ \\ z[ - v. 4^^^^' ' 

(5.75) 

where the contour Fp encloses = ^^/(l — qi) anticlockwise, and Cr^ encloses the origin 
z' and Pi = gj/(l — qi). Introducing yiji) {% = 1, 2), which will be chosen later, we set 

k(ti,Xi;t2,X2) 

^ / ^ / "^2 ^2 T(hr,Az'o)-hr, jz'-, )) TT ^2 ~ -Pt ^1 / + ^1 \ 

i2my ^ 4 4 4 - 4 11 4 - ,,^^^-yit.)T+M yi + z'J ' 

(5.76) 

where the parameter r is defined in (I3.54p and 

hr{z) = e^" ln(l + 4 - y{T) \n{z). (5.77) 

We choose ?/(r) in a way that hT-{z) has the saddle point at z = p = q/{l — q) where q is 
defined in (I3.56p . Thus from the condition, h'^{p) = 0, we have 

y{r) = (5.78) 
l+p 

We scale the variable 4 = 1) 2) around the saddle point Zc = p as 

..4l + ^l(i±P)J. (5.79) 
V en(pT)' ) 

Considering (13.550 . we also scale Xi (i = 1, 2) as 

Xi = ——T + ii^ • (5.80) 

From these equations and (13.560 . we find 

' 4 - 4 6-^2^2 - e-^^wi V2(l + p) ' 



rxi-y{Ti)T 



^1 + 4^*''-^ 
1 + 4 



4 - - - ^1 - Pi jL- e -u/i -r tj 

1. (5.81) 
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From these equations, we eventually obtain 

(5.82) 

Here the contour F encloses — e'^^ej, (i = 1, ■ ■ ■ , M) anticlockwise and 7 is arbitrary path 
from —ioo to ioo. 

Next we consider the scaling limit of 4'ti,t2{xi,X2) (13. 7p . Changing the variable z to 

2;' = 1/2;, we set 

^t,Mi^i,X2) = — /^e^('^-^(^')-'^^i(^'«^^-^^-(^(-^)-^^(^^))^, (5.83) 
27rz / z' 

where hr{z) and ?/(r) are given in (15.771) and (15.781) . We set 



From this equation and (15.801) 

^T{hr2{z)-hr-,{z)) ^ ^T{hr^{z)-hr^{z)) ^W^ 



^„,-iy^,-y^,)T ^ / ^ (^^ " ^^^^ ) . (5.85) 



Combining these equations, we finally find 

0ti,t2(a;i,X2) 



,t.-t.^.-.. v^(i +p) 1 r .... „„„ (...2 , 2(6 - eo^ 



- ^^+^)^^-^^^^^"^'^^V .(l-e2(-^-n)) -P I T-e2(n-.) J • (^-SS) 

Thus from (15.821) and (15.861) and noting that the prefactor (1 +p)*2"*^p^^~^2 does not affect 
the determinant, we get the desired expression for the kernel. 



6 Discussion 

6.1 Numerical simulations 

In this section, we give the result of Monte-Carlo simulations about the position fluctuation 
of a tagged particle in the TASEP in order to check our analysis discussed in the preceding 
sections. 
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We performed the simulations in two situations, the case without defect particle and 
that with one defect particle. Fig. 8 shows the data of the scaled position of a tagged 
particle obtained by the Monte-Carlo simulation and the probability distribution functions 
which must fit them. Fig. 8(a) corresponds to the time region 1/(1 — q) < u < Uc = 
{(f — 2qq + q)/{q — qY where q (resp. q) is the stay rate of the normal particle (resp. 
the defect particle). Note that in the region, the both situations belong to the region 2 
and due to Theorem 2-2, the GUE Tracy- Widom distribution must fit the data. Fig. 8(b) 
represents the case u = Uc- In this case the first situation comes under the region 2. On 
the other hand, the second one belongs to the region 3 where the position fluctuation are 
described by the limiting largest eigenvalue distribution of GOE^ as explained in Theorem 
2-3. In Fig. 8(c), only the data for the second situation are shown for the case Uc < u. This 
case is classified as the region 4 and from Theorem 2-4, the fluctuation is supposed to be 
Gaussian. In all figures, we see a good agreement between the data and the distribution 
functions which they must obey. 



6.2 Correlation of current fluctuations 

In this article, we have discussed the multi-time distribution of the tagged particle (12.11) . 
In order to analyze it, we have introduced the directed polymer problem of the 01 matrix 
in Section m 

Here we consider the directed polymer problem of another random matrix where each 
element is geometric random variable. By this analysis we can discuss the correlations of 
other quantities in the TASEP with the step initial condition. 

Let {ac,(z, j)}i<j<Ar,i<j<M be the N x M matrix and the element ag{i,j) is the geometric 
random variable, 

Prob(a,(z,j) = A;) = (l-g,)gJ, (6.1) 

where qj is the parameter of the geometric distribution and we identified this with the stay 
rate of the jth particle of the TASEP with the step initial condition defined in Section [2l 
For the matrix, we introduce the quantity G*g{N,M) analogous to G{N,M) (14.31) as 

G*g{N,M) = N + M -1+ max ^ ag{i,3). (6.2) 
Here T^g{N, M) is the set of right /down paths from (1, 1) to (A^, M), 

■Kg{N,M) = {{{ik, jk)k=l,2,- ,N+M-l}\{il, jl) = (1, 1), (^AT+M-l, jAT+M-l) = {N,M), 

{ik+1 -ik,jk+i -jk) = (1,0) or (0,1)}. (6.3) 
In this setting, we consider the following quantity, 

Prob {G*g{N,, M) < ti, • ■ ■ , G*g{N^, M) < Q . (6.4) 
It is known that G*g{N, M) can be interpreted as a quantity in the TASEP as follows pTl 
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G*g{N,M) = t: In the TASEP with the step initial condition, the time until 
which the Mth particle moves N sites to its right is t. 



Thus the quantity (16.41) represents the correlations of the arrival times of the Mth particle 
at the site Ni {i = 1, ■ ■ ■ yfn). (Note that we set the site coordinate as in Fig. 1(b).) 
Furthermore we easily find that 

, /^*. , , ( Fioh (Hit, N) > N) , foriV<M, , , 

Prob G* A^, M < t) = <^ ~ ' ~ 6.5 

^ ^~ ^ [Pioh {H{t, N) = M) , for M<iV. ^ ^ 

Here H{t, N) represents the number of particles which passed the site N until time t. Thus 
the probability (16.41) also means the correlations of currents between different times and 
different sites. 

Applying the similar technique in Section H] to (16.41) . we can also represent it as the 
growth process of Young diagram characterized by the Schur process. The result is as 
follows, 

Prob {GliN,, M) < ti, ■ ■ ■ , GliN^, M) < 

= Prob (XiiNi, M) <ti- N - M + !,■■■ , Ai(A'„, M) < t„, - N - M + 1) , (6.6) 

where A(A^j, M) is the Young diagram obtained by applying the normal RSK algorithm to 
the submatrix {ag{i, j)}i=i^... ^N„j=i,--- ,m and Ai(A^i,M) is the length of its first row. The 
probability measure in the right hand side of this equation is characterized by the following 
joint distribution function. 



M 

N 



Prob (A(l, M) = A(^), ■ ■ ■ , XiN, M) = A^^^) 

= Sx(N){qi,q2, ■ ■ ■ ,gM)sA(i)(l,0, ■ ■ ■ )sa(2)/a(i) (1, 0, ■ ■ ■) " ■ ■ Sx(n)/x(n-i){1,0, ■ ■ ■) " 

1=1 

(6.7) 

where sx/^{qi, q2,- ■ ■) is the Schur function. This is also a special case of the Schur process 
and we can get the Fredholm determinant representation of (16. 6p . Although in one point 
case, the asymptotics of this equation was discussed in [IZlllS], the multi-point distribution 
has not been discussed yet. 



7 Conclusion 

In this article, we have studied the multi-time distribution function (12. ip of position fluctu- 
ations of a tagged particle in the TASEP with the step initial condition. The main results 
are summarized as Theorems 1, 2 and 3 in Section [3l 

First, we have obtained the Fredholm determinant expression of (12. ip in Theorem 1. 
For this purpose, we have mapped the time evolution of a tagged particle in the TASEP 
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to the growth process of Young diagram which is related to the special case of the Schur 
process. Next, using the Fredholm determinant in Theorem 1, we have studied the two 
types of scaling limit. The first one is the case where both time t and the label of a tagged 
particle M go to infinity. In the second one, we take the t ^ oo limit with M fixed. 
The results for the first and second ones are shown in Theorems 2 and 3 respectively in 
Section [31 

In the first scaling limit, if the hopping rates of all particles are the same , we can divide 
the scaled time into two characteristic regions according to the limiting behavior of (12.11) . 
the region 1 where a tagged particle begins to move and the region 2 which is after the 
region 1. In the region 1, the limiting process of a tagged particle converges to the spatially 
discrete process which is described in Theorem 2-1. The process reflects on the discreteness 
of the model. In the region 2, on the other hand, we have shown in Theorem 2-2 that the 
limiting process becomes the Airy process, which is characteristic of the one-dimensional 
KPZ universality class. 

If there are n defect particles with small hoping rates in front of a tagged particle, the 
limiting distribution changes at the scaled time Uc which is determined by the hopping 
rate of the slowest defect particle. The limiting process around Uc (region 3) is equivalent 
to that of the largest eigenvalue in GUE Dyson's Brownian motion with rank one external 
source as described in Theorem 2-3. Theorem 2-3' indicates that the rank becomes n if the 
hopping rates of the defect particles are the same. When the scaled time is after Uc (region 
4), we have found in Theorem 2-4 that the process is equivalent to the one dimensional 
Brownian motion. If the hopping rates are degenerate, it is equivalent to the process of the 
largest eigenvalue of n x n GUE Dyson's Brownian motion model as described in Theorem 
2-4'. This indicates that in the region 4, the effect of the defect particles is dominant 
whereas that of an infinite number of normal particles is irrelevant. 

Theorem 3 shows the result for the second scaling limit. The limiting distribution is 
also described as the Fredholm determinant whose kernel is the same as that in Theorem 
2-4'. 
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Figure Captions 



Fig. 1: Definition of tlie TASEP witli tlie step initial condition, (a) During eacli time step, 
tfie ith particle can hop to the right with probability 1 — g^. However, i + 2th particle 
cannot hop to the right neighboring site since the site is occupied by i + 1th particle, (b) 
Step initial condition. 

Fig. 2: Typical time evolution of the particles from the first to 100th in the TASEP with 
the step initial condition. The hopping rate of all particles is 0.9 (g = 0.1) except the 
four defect particles (the first, 25th, 50th and 75th particles) whose hopping rate is 0.8 
(g = 0.2). Fig. (a) illustrates the whole time evolution from t = through 3000. Its closeup 
around t = 200 and t = 3000 is shown in (b) and (c) respectively. 

Fig. 3: Average position of the tagged particle. The thick line shows (13.111) . We divide 
this line into four regions, where u is near 1/(1 — q), 1/(1 — g) < u < Uc, u is near Uc, and 
Uc < u. We denote them as region 1, 2, 3, and 4 respectively. 

Fig. 4: Example of the time evolution of the TASEP and corresponding 01 matrices. In this 
case, two tables is assigned for one example of the time evolution. Since the second particle 
can not hop to its right neighboring site between t = 2 and 3, there are two possibilities of 
assigning both and 1 in the (i, j) = (1, 2) element of the 01 matrix. 

Fig. 5: (a) Example of time evolution of the TASEP and (b) one of its corresponding 
matrices. Note that a set of many matrices corresponds to (a) and only one of them is 
shown here. 

Fig. 6: (a) Generalized permutation corresponding to Fig. 5(b). We arrange where 
o-ihi) = 1 following the rule in [60]. (b) The pair {P,Q) obtained from (a) by the dual 
RSK algorithm. 

Fig. 7: Growth of Young diagram obtained from Fig. 5(b) by the dual RSK algorithm. The 
diagrams A(i,4) (1 = I,-- - ,6) are obtained from the submatrices {a{j, k)}j=i^... ^i^k=i,--- ,4 
of Fig. 5(b). The set of the diagrams {A(i, 4)}j=i_... corresponds to the Q tableaux in 
Fig. 6(b). 
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Fig. 8: Probability distributions of the scaled position of the 100th particle from the right 

(M = 100) for t = 200 (a), t = 1000 (b) and t = 3000 (c). In these figures, x represents 
the data for the case without defect particles and we set the stay rate q = 0.1. On the 
other hand, + corresponds to the case where the first particle is a defect particle with 
q — 0.2 while remaining particles are normal ones with q = 0.1. The number of samples 
are 10000 for each case. In (a), both cases belong to the region 2 and are fitted into the 
GUE Tracy- Widom distribution shown as the dashed line. In (b), the second case belongs 
to the region 3 since u = t/M = = u^. where u^. = {(f — 2qq + q)/{q — qY ^-nd it is 
described by the distribution denoted as GOE^ while the first case remains in region 2. 
In (c), the second cases comes under the region 4 where the distribution is described by 
Gaussian (dashed line). 
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Fig. 6 
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